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Abstract 

The long-time relaxation of ideal two dimensional (2D) magnetohydrodynamic 
^ , (MHD) turbulence subject to the conservation of two infinite families of constants of 

CN I motion — the magnetic and the "cross" topology invariants — is examined. The analysis 

■ of the Gibbs ensemble, where all integrals of motion are respected, predicts the initial 

. state to evolve into an equilibrium, stable coherent structure (the most probable state) 

I and decaying Gaussian turbulence (fluctuations) with a vanishing, but always positive 

temperature. The non-dissipative turbulence decay is accompanied by decrease in both 
the amplitude and the length scale of the fluctuations, so that the fluctuation energy 
remains finite. The coherent structure represents a set of singular magnetic islands 
with plasma flow whose magnetic topology is identical to that of the initial state, while 
the energy and the cross topology invariants are shared between the coherent structure 
, and the Gaussian turbulence. These conservation laws suggest the variational principle 

of iso-topological relaxation which allows us to predict the appearance of the final state 
from a given initial state. For a generic initial condition having x points in the mag- 
netic field, the coherent structure has universal types of singularities: current sheets 
^ ' terminating at Y points. These structures, which are similar to those resulting from 

d • the 2D relaxation of magnetic field frozen into an ideally conducting viscous fluid, are 

observed in the numerical experiment of D. Biskamp and H. Welter [Phys. Fluids B 
1, 1964 (1989)] and are likely to form during the nonlinear stage of the kink tearing 
mode in tokamaks. The Gibbs ensemble method developed in this work admits ex- 
tension to other Hamiltonian systems with invariants not higher than quadratic in the 
highest-order-derivative variables. The turbulence in two dimensional Euler fluid is of 
a different nature: there the coherent structures are also formed, but the fluctuations 
about these structures are non-Gaussian. 
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I Introduction 



Turbulence in two dimensional fluids is simple enough to allow some analytically tractable 
models [|l], ^, ||, ^ and high-resolution computation [H, ||, 0, ||, , but still sufficiently complex 
to exhibit the general challenging nature of turbulence. In addition, many features of plasma 
and geophysical fluid dynamics are essentially two dimensional. 

The principal difficulty encountered by any analytical approach to fluid turbulence is 
that the underlying equations of motion are nonintegrable, a quality independent of the 
abilities of the researcher. To overcome this difficulty, two major analjd;ical approaches have 
been attempted so far, namely, (a) statistical closures dealing with simplifying, but usually 
uncontrollable, modifications of the dynamical equations and (b) considerations of a priori 
equally probable states constrained by the known integrals of motion. In this work we take 
the second approach, where the consistency requirement is to allow for all constants of the 
motion. 

The conventional Boltzmann-Gibbs statistical mechanics, when applied to partial differ- 



ential equations |10, 11, 12, 13], encounters a fundamental obstacle because of the underlying 
infinite number of degrees of freedom and the necessity to use a finite (A^) dimensional ap- 
proximation. The continuum limit N —>■ oo, in common to all classical fields, results in the 
so-called ultraviolet catastrophe (the divergence of energy at finite temperature), a problem 
which goes back to Jeans. Unlike the equilibrium electromagnetic radiation, the ultraviolet 
catastrophe in fluid turbulence cannot be remedied by quantization and should be resolved 
within the classical framework. The tendency toward the equipartition of energy between 
the degrees of freedom (by T/2 for each degree, where T is the energy temperature) in a 
closed continuum system can only be satisfied by letting the temperature to zero. In fact, 
the way the temperature goes to zero as the number of degrees of freedom goes to infinity 
is the heart of the problem. Nontrivial equilibrium states are obtained when there are more 
than one integrals of motion, which diverge at different rates a.s N oo. 

The absence of a well-defined concept of measure in a functional (infinite dimensional) 
space requires an dimensional discretization, even though the final results are obtained 



by letting — oo. Lee [Tj] was the ffist to use the truncated Fourier series to show the 
validity of an infinite dimensional Liouville theorem. The choice of the discrete variables is 
not unique, and one ought to make sure that the results of the statistical theory of turbulence 
be invariant with respect to the way this choice is made. 

A less formal, physical motivation for the discretization procedure can be found in the 
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finiteness of the number N{t) of the "effectively excited" collective modes at any finite time 
t. Under "effectively excited" we mean modes with amplitudes not yet exponentially small 
(a smooth field, for instance, has exponentially small high Fourier harmonics). There are 
examples of N{t) becoming infinite in finite time [1^, even when a real-space collapse |jl6 



does not occur. It appears, however, that in two dimensional ideal hydrodynamics and 
magnetohydrodynamics N{t) behaves algebraically, N{t) ~ {t/TAY, so that the time of 
the doubling of N{t) is of order t and the time of the increment of N{t) by one is of order 
6t = TA{t / ta)^~^ , if < p < 1. Hence, on time scale longer than the eddy turnover (nonlinear 
mode interaction) time r^, one may expect an equilibrium statistical distribution among the 
effectively excited N{t) modes. By letting N{t) ^ 1 in this distribution, we shall infer the 
most probable direction of the long-time system evolution. The probability of a deviation 
from the most probable state goes to zero as the number of degrees of freedom goes to 
infinity. As for a continuum conservative system we really mean N{t) ^ oo as t oo, the 
probabilistic nature of the statistical prediction assumes a rather deterministic quality. 

Artificial finite dimensional approximations are notorious for destroying an infinity of 
topological invariants (also known as freezing-in integrals or Casimirs), which impose im- 
portant constraints on the evolution. An interesting alternative proposed by Zeitlin ||17|| , 
whereby an N dimensional hydro dynamic-type system conserves ~ invariants, is not 
quite suitable for continuum fields because of the implied periodicity in the Fourier space 
corresponding to modulated point vortices in real space. It would be very interesting to 
construct other "meaningful" (that is, having many invariants) finite-mode hydrodynamics 
with well-behaved real-space velocity fields. 

So far, most statistical theories of continuum hydrodynamics, most notably the absolute 
equilibrium ensemble (AEE) theory 0, ^], simply ignored all topological invariants 
except quadratic ones, such as enstrophy and helicity in hydrodynamics or magnetic helicity, 
cross helicity, and square vector potential (in two dimensions) in magnetohydrodynamic 
turbulence. These integrals were honored the special attention in part because of their 
ruggedness (survivability under the approximation of a finite number of Fourier modes), but 
mostly because of the convenience to handle quadratic integrals. 

Despite the dissatisfaction with such a reasoning (cf. [0), the attempts to incorporate 
all topological invariants in Gibbs statistics have been less frequent, the examples including 
Vlasov-Poisson system and 2D Euler turbulence p, |2^, |2^, |2^. For the reasons dis- 
cussed in Sec. |VJ and Appendix D, these attempts appear not quite successful, because the 
non-Gaussianity of turbulent fluctuations in these systems poses a fundamental difficulty in 
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making quantitative predictions about the equilibrium turbulent states. The key problem is 
that the averaging with respect to the given probability functional of a turbulence involves in- 
tegration in a functional (infinite dimensional) space. These averages are well defined — that 
is, independent of the discretization procedure involved, if and only if the probability func- 
tional is Wienerian, or Gaussian in the highest-order derivative [^. Otherwise, the result of 
the averaging is sensitive to the arbitrary choice of the sequence of discrete representations, 
which makes such probability functionals ambiguous, to say the least. In the present work 
we point out that these difficulties are absent from a certain class of systems; namely, those 
where all integrals of motion are not higher than quadratic in the highest-order-derivative 
variable. Important examples of such systems, allowing a valid Gibbs-ensemble description 
of turbulence, include two dimensional and reduced magnetohydrodynamics. 

The problem of accounting for all invariants is circumvented by the representation of 
turbulence in the form of a gas of point vortices, which is a very singular, and very special 
topologically, although an asymptotically exact solution of the hydrodynamic equations. The 
localization of vorticity in point vortices makes the topological constraints trivially fulfilled 
for any motions. The conservation of only energy and the Liouville theorem expressed in 
the convenient form of the spatial variables of the vortices yield nicely to the statistical 
mechanical description, although the thermodynamic limit of infinitely many point vortices 
has long been a controversial issue [0, ^ |2^, |3^. It remains unclear to what extent 
the gas of many point vortices represents a continuum two dimensional turbulence. The 
frustrating dependence of the statistics of point vortices on the arbitrary choice of their 
strengths was noted by Onsager and refiects the above-mentioned fundamental difficulty in 
the 2D Euler turbulence. 

We pursue an analytical approach to continuum two dimensional ideal MHD turbulence 
where all topological invariants are respected. We use the Gibbs ensemble analysis to predict 
the following evolution of the turbulence. An initial state evolves into (a) a stationary, stable 
coherent structure, which appears as the "most probable state" and (b) small-scale turbu- 
lence (fiuctuations) with Gaussian statistics. The Gaussianity of the MHD fiuctuations was 
recently numerically confirmed by Biskamp and Bremer |3l|. At large time the fiuctuations 
of the magnetic fiux and the fiuid stream function assume vanishing amplitude and length 
scale (while containing finite energy and dominating phase volume), and become essentially 
invisible on the background of the coherent structure. In this sense, the coherent structure 
can be regarded as an attractor or a relaxed state, although the underlying dynamics is 
perfectly Hamiltonian. 
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The concept of "statistical attractor," introduced by Vladimir Yankov p2|, |33|, 



em- 



phasizes the method of analysis and describes this kind of Hamiltonian relaxation, when the 
excess of phase-space volume and energy get hidden in obscure (small-scale) corners of the 
infinite dimensional phase space. The fundamental difference between the statistical attrac- 
tors in nonintegrable wave-type systems |33|, |3^, ^ (where there are only a finite number 
of integrals of motion) and the hydro dynamic-type systems (where the number of integrals 
is infinite) is the universal shape of the attractor — soliton — in the first case, and the nonuni- 
versal shape of the coherent structure — vortex — in the second case. The appearance of the 
coherent structure depends on the initial condition, but is the same for all initial conditions 
with identical topological invariants [Eqs. (6) and (7)]. In this sense, the asymptotically 
emerging coherent structures (relaxed states) in hydrodynamic systems are attractors only 
within a certain subclass (basin) of initial conditions. These relaxed states can be called 
topological attractors. 

The specific of two dimensional MHD is that the coherent structure inherits from the 
initial state all magnetic topology invariants, but only a fraction of energy, the rest of which 
goes to the Gaussian fluctuations. These invariant sharing properties can be interpreted in 
terms of the well-known reasoning of turbulent cascades. In the case of 2D MHD the energy 
cascade is direct (i.e., toward the small-scale fluctuations), while the magnetic topology 
cascades inversely (toward a large-scale magnetic structure). However, compared to the 
cascade description, the invariant sharing properties appear to present a clearer physics 
of what happens to the conserved quantities in a closed turbulent system. In fact, this 
allows us to predict the appearance of the relaxed state, which should minimize energy 
subject to certain topological constraints. One of the novelties of our analysis is using a 



second functional set of "cross" topology invariants [^, |3^, which was not used in MHD 
turbulence theories so far, including the previous note |3^. We find that these invariants 
have an important effect on the statistical description of turbulence; specifically, the shape 
of the coherent structure is sensitive to the cross topology invariants. In many features 
our development is analogous to the Gibbs ensemble treatment of the truncated Fourier 
representation of 2D MHD which was investigated by Fyfe and Montgomery pO[. The reason 
for this resemblance is partial consistency of the Gibbs statistics with partial invariants 
accounting. In fact, our theory shows that truncated 2D MHD equations partially represent 
the true statistics of ideal 2D MHD. 

Unlike neutral fluid dynamics, magnetohydrodynamics in two dimensions are known to 
produce energetic small scales. This makes the difference between 2D and 3D for MHD 
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turbulence less drastic than that for fluid turbulence. In our model, we observe two types of 
small-scale behavior: (a) for a generic (that is, topologically nontrivial) initial condition, the 
coherent structure must have discontinuities in the form of current sheets and (b) the Gaus- 
sian fluctuations in the long evolved state have both vanishing length scale and amplitude 
so that the gradients and the energy are flnite. The numerically observed current- sheet-type 
structures are explained by our theory as the singular coherent structures. 

The paper is organized as follows. In Sec. |I| we present the equations and the constants 
of the motion. In Sec. ^ we discuss the stationary solutions of the MHD equations and 
formulate the Arnold variational principle in the form suitable for 2D MHD. In Sec. ^ 
the canonical ensemble approach to MHD turbulence is set forth (Sec. [IV. Ap , and the most 
probable state (Sec. |1V.BD and the fluctuations about this state (Sec. |1V.(JD are analyzed. 
The key issue of how the integrals of motion are divided between the coherent structure 
and the fluctuations is addressed in Sec. |IV.D| . In Sec. we reformulate the properties 
of the coherent structure using a variational principle of iso-topological relaxation, which 
allows us to predict the appearance of the structure. This prediction is then compared with 



numerical results (Sec. [V.A|) . In Sec. [V.B| we speculate on the role of small dissipation and 
the relation between the resistive and the ideal MHD relaxation in the kink tearing mode in 



tokamaks. Section ^ restates the principal steps of our statistical method and summarizes 
our work. Some technical details and results not directly related to MHD turbulence are 
set in Appendices. In Appendix A we discuss the Lyapunov stability of MHD and Euler 
fluid equilibria and point out the relation of minimum- and maximum-energy stability to 
positive- and negative-temperature Gibbs states, respectively. Appendix B addresses the 
Liouvillianity of the eigenmodes that we use in the Gibbs statistics. In Appendix C the 
spectrum of the eigenmodes is studied. In Appendix D we discuss the application of the 
Gibbs-ensemble formalism to the turbulence of two dimensional Euler fluid. 



II Equations and constants of the motion 

We consider the set of equations of two dimensional incompressible ideal magnetohydrody- 
namics (cf. [09| ) 

dta = {il},a} , (1) 
dtu = {^,uj} + {j,a}, (2) 

where {A, B} = x V-B-z denotes the Poisson bracket, z the unit vector in the z direction, 
ip{x,y,t) the stream function of the fluid velocity fleld v = V x (ipz), u = —'V'^ip the fluid 
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vorticity, a{x, y, t) the normalized vector potential of the magnetic field B = (47rp)^/^ V x (az) 
(p being the constant fiuid density), and j = — V^a the normalized electric current fiowing 
perpendicular to the {x, y) plane. The boundary conditions a = ip = ^ are assumed at a 
rigid boundary encompassing the finite domain with the area S. 

The incompressibility of the fiuid motion is a reasonable approximation in tokamaks 
where the strong toroidal magnetic field makes plasma compression energetically expen- 
sive. In the reduced MHD approximation where the fast Alfven and magnetosonic 



waves due to are ignored, the strong uniform field B^ drops out of the equations of the 
motion. 

The system (1) — (2) conserves the following quantities: the energy 

E = E^ + Ef = ^J [(Va)^ + (V^py] d^^=lj («J + V^^) ^^'x , (3) 

consisting of the magnetic part and the fiuid part Ef, the momentum (with translation- 
ally invariant, or in the absence of, boundaries) 

V ci^x = - f XX zuj (fx , (4) 



2 

the angular momentum (with circular or no boundaries) 



Mz = / x X vfi^x = - / X X (x X zu) £x , (5) 



3 

the magnetic topology invariants 

Ip = j F{a) (fx , (6) 

and the "cross" topology invariants 

Jg = j ^G{a) (fx = J G'(a) • Va (fx , (7) 

where F and G are arbitrary functions. Along with the continuum set of integrals (6) and 
(7), we will also use their discretized analogues. 

In = J a"" (fx , (6a) 
J„ = y ua!^ dfx , (7a) 

through which the continuum invariants can be Taylor expanded. 
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Strictly speaking, invariants (6) and (7) do not yet imply the conservation of topology. 
Equation (6) only means that allowed motions are incompressible interchanges of fluid el- 
ements together with their "frozen" values of the magnetic flux a. The topology of the 
contours of a will be conserved only if these interchanges are performed by continuous move- 
ments, a constraint which is not built into Eq. (6) but follows from the equations of motion 
for smooth initial conditions. Then the conservation of magnetic topology expressed by 
Eq. (6) means that if the contour a{x,y) = ai initially lies inside the contour a{x,y) = 02, 
then this topological relation is preserved by the motion; if a contour of a has a hyperbolic 
(saddle, x) point, this quality will also persist. In addition to the topological constraints, 
the integrals (6) also specify the incompressibility of the fluid, so that the area inside a given 
contour of a remains constant. The geometrical meaning of the integrals (7) is that the 
amount of the fluid vorticity u on a given contour of a (to be more precise, the integral of u 
over the anulus between two infinitesimally close contours) is conserved. Although stating 
nothing of the contours of 00 or tp, the conservation of the integrals (7) also bears certain 
topological relation between the magnetic field and the vorticity field, which motivates our 
notation of the "cross topology invariants." 

The invariants Jq appear to be poorly known, although the particular member of the 
family (7a) — the cross helicity 



Ji = [ ua (i^x = —== / V ■ B d^x 



— has been extensively discussed in the literature. The invariants (7) were first noted by 
Morrison and Hazeltine |3^. Independently, a similar set of integrals was used to study the 
vortex stability in the framework of two dimensional electron MHD 1^^. The idea towards 
the existence of a second set of topological invariants is suggested by the observation that 
there is another frozen-in quantity, namely the vorticity u, in the Euler limit a = 0, which 
must have a counterpart in the MHD case a 7^ 0. Once the existence of a second functional 
set of invariants is suspected, it is not hard to guess the form of the topological invariants (7). 
Although this is not straightforward, one can trace the transition, as a — 0, from Eqs. (6) 
and (7) to the Euler invariants / F{uj) (i^x. 

Another way to find the topological invariants is to identify the Hamiltonian structure 
using a noncanonical Poisson bracket ||36|, whereby the topological invariants appear as 
Casimirs. 
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Ill MHD equilibria and Arnold's variational principle 



The system of equations (1) — (2) has stationary solutions satisfying 

{^,a} = 0, (9) 
= {a,VM- (10) 

The equihbrium condition can be rewritten in a more convenient form by substituting the 
functional dependence = ipla), which is implied by Eq. (9), into Eq. (10). Then, after 
simple manipulations, we find 

{^{a),V'^{a)} = , (11) 

where 

^(a) = £ da'^± (1 - [d'il){a')/da'Y) , (12) 

and the sign is chosen to make the square root real. As equation (11) shows, any two 
dimensional MHD equilibrium with fluid flow {ip ^ 0) is reduced to a purely magnetic 
[if) = 0) equilibrium for a modified magnetic vector potential a' = ^^(a). Note that the 
magnetic field lines (the contour lines of the vector potential) are identical for both the true 
field a and the modified field a', although the values of a and a' on these lines are different. 

The most evident stationary solution is given by an arbitrary circular distribution ip = 
ip{r), a = a{r), where r is the distance from the origin. Another solution to (9) — (10) 
corresponds to the identical zero in one of the Elsasser variables, where il){x,y) = ±a{x,y) 
is an arbitrary function of x and y, so that "^{a) = 0. For the case when {tpl and \a\ are not 
identical, there exist many periodic and quasiperiodic solutions in the form 

N 

C^(x) = a(x) = J2 cos(k^x + Om) , (13) 

m=l 

where the moduli of the wavevectors are the same. In addition to the smooth solutions 
one can devise a wide class of singular solutions with appropriate boundary conditions at the 
lines of discontinuity. Without additional physical constraints, such as stability or topology, 
the class of all equilibria is too wide to be useful. In Sees. |3 and we provide such 
constraints, which specify the physically interesting (attracting) equilibria. In many cases 
these equilibria must be singular. 

There exists a profound relation between the stationary solutions and the constants 
of the motion. For finite dimensional conservative systems, the D'Alembert variational 
principle says that the energy variation be zero at an equilibrium. For a hydrodynamic-type 
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conservative system, the counterpart of the D'Alembert theorem is the Arnold variational 
principle. Originally formulated for the Euler equation, but carried over without difficulty to 
other hydrodynamics, the principle states that at a stationary (and only stationary) solution 
the variation of energy, subject to the conservation of all topological invariants, be zero: 



SE 



= . (14) 

Ip ,jQ=const, V-F and G 



If, in addition, the second variation is definite — that is, the energy assumes a nondegenerate 
conditional extremum, then the equilibrium is Lyapunov stable. In fact, this was the search 
for stable fiuid fiows which motivated Arnold's work. The power of the method lies in the 
possibility to write the general iso-topological (iso-vortical in Arnold's notation) variation in 
a closed form. When all the integrals (6) and (7) are to be conserved, such a variation is 

where the infinitesimal iso-topological variation is given by |37 



/i(x, y) and v{x,y) being arbitrary functions. The operator of finite variation A can be 
symbolically expressed through the infinitesimal variation 5 as 

A = exp5-1. (17) 

It is easy to verify that the finite variation (15) conserves both sets of integrals Ip and Jg to all 
orders in /i and v. The form of the variation is suggested by the form of Eqs. (1) — (2), where 
one can substitute the quantities if) and j by arbitrary dfi/dt and du/dt, respectively, to 
preserve only the Poisson-bracket structure of the equations and thereby to iso-topologically 
(that is, at constant Ip and Jg) drag the fields a and t^; to a new state, where the values of 
all other integrals, if any, are generally different from those of the initial state. 

Let us see what happens to the energy (3) under the variation (15) — (16). Writing the 
total change in the energy in the form 6E + 6'^E/2\ + . . ., we obtain after integrating by parts 

SE=-J[fi ({^, uj} + {j, a}) + iy{^, a}] rf^x . (18) 

By requiring that the first energy variation (18) be zero for all n and u we arrive exactly at 
the system of equations (9) — (10) specifying the equilibrium solution. This strongly suggests 
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that the two sets of the topological invariants (6) and (7) are indeed complete, a condition 
necessary to apply statistics (Sec. fV]) in a meaningful way. 

The second iso-topological variation of energy can be used to investigate the stability of 
MHD equilibria, as discussed in Appendix A. 

IV Gibbs statistics of two dimensional MHD turbu- 
lence 

We now wish to analyze the long-time evolution of the system subject to Eqs. (1) and (2) for 
a given initial condition ao(x), ipo{^)- The probability distribution functional P[a(x), ^/^(x)] 
can serve this purpose. This functional specifies the relative probability, with respect to 
time measure, of the spatial behaviors of various states a(x, t) and '?/'(x, t). As P is invariant 
under the evolutional change of the fields a and tp, it must be a function of the constants of 
motion (3) — (7). 

IV. A Choice of statistical ensemble 

For a conservative Hamiltonian system, P is given by the microcanonical ensemble, 

PMc[a,^] = 5{E[a,^P] - E^) H^UnK^] - /„o) 5(^n[a, ^] - J„o) , (19) 

n 

specifying a uniform distribution on the manifold of the specified (initial) integrals of motion. 
It must be emphasized that the validity of the microcanonical ensemble requires at least four 
assumptions. 

1. The phase space must be finite dimensional; that is, the fields a and ip are parameterized 
by discrete dynamical variables /m, m = 1, 2, . . . , A^, so that the concept of measure 
in the space of states is meaningful. This is a tricky issue as to how many variables 
are needed (see discussion in Sec. |). 

2. The motion on the manifold of conserved invariants is nonintegrable (chaotic). The 
fundamental phenomenon lying behind the ergodic behavior (uniform distribution over 
the manifold) is Hamiltonian chaos, whose principal manifestation is the exponential 
divergence of nearby trajectories. The chaotic motion is possible only if the dimension 
of the manifold (A^ — A^/, where Nj is the number of invariants) is three or more (four 
or more to allow the Arnold diffusion, so that all of the manifold might be visited 
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by each trajectory). The property of ergodicity was proved for very special cases |^ 
but is beheved to be generally valid if the dimension of the manifold is sufficiently 
large: N — Nj ^ 1. The remarkable accuracy of the classical thermodynamics is 
connected with the macroscopic numbers of degrees of freedom (A^ ~ 10^^) and only a 
few invariants. It is also natural to expect that the microcanonical statistics will work 
in turbulence (which can be loosely defined as chaos in PDE), where the limit N oo 
should be carefully taken. 

3. The manifold of conserved integrals specified by the delta functions in (19) must be 
connected. The problem of connectivity of complicated iso-surfaces is related to the 
percolation problem |Q. The percolation threshold, above which the connectivity 
takes place, is inversely proportional to the dimension of the iso-surface. This suggests 
that in the continuum (infinite dimensional) limit the connectivity of the manifold of 
conserved integrals should not be a problem. 

4. The dynamical variables fm(t) must satisfy the Liouville theorem: J2mdfm/dfm = 0. 
For non-Liouvillian variables a weighting factor (the Jacobian of change to Liouvillian 
variables) should be included in Eq. (19). 

In a hydrodynamic-type system, where the number of dynamical constraints is infinite, 
we encounter another difficulty. Namely, any attempt to restrict the dimension of the phase 
space without restriction on the number of conserved quantities immediately drives the 
manifold of conserved integrals into an empty set, where no mixing may occur. Motivated 
by the experimental/numerical observation that turbulence does exist, as well as by the 
functional arbitrariness of the iso-topological variation (16), we adopt a hypothesis that there 
exists a "meaningful" N dimensional MHD approximation with at most Nj{N) conserved 
integrals, where both Nj and N — Nj go to infinity as ^ oo. (In Zeitlin's example ]l^ 
for the Euler fiuid Nj ^ N^^"^ .) The specific form of this approximation is unimportant for 
our arguments. 

The microcanonical ensemble (19) is inconvenient to handle and is commonly transformed 
into the more convenient canonical (Gibbs) ensemble by integrating Pmc over most of the 
dynamical variables in the amount of Nth ^ N — Nth ^ 1- These Nth degrees of freedom 
can referred to as "thermal bath." The integration over the thermal bath variables leads to 
an exponential dependence of the resulting distribution on the integrals of motion expressed 
through the remaining N — Nth variables, the rest of the information being stored in arbitrary 
constants called temperatures. 
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In our problem, the dimension N' = N — Nth of the subsystem can be also taken large, 
which amounts to another finite dimensional (A^') MHD approximation. However, now we 
have the canonical distribution over the remaining A^' variables. 



P[a, ip] = exp 



(20) 



instead of the microcanonical one (19). A drastic simplification achieved by the change of 
the ensemble is that in the finite dimensional Gibbs distribution (20), where the fields a and 
ip are parameterized by A^' modes and the summation over invariants runs up to Nj{N'), 
we may extend the summation up to infinity without significant change in the result, which 
was impossible for the product of delta functions in the microcanonical ensemble (19). 

In Eq. (20), the constants a, a(3n and 07^ appear as the reciprocal temperatures corre- 
sponding to each invariant. These constants are to be determined from the initial state by 
solving the infinite system of equations: 

{E)p = Eo, (/„)p = /„o, {Jn)p = Jno, n = 0,1,2,... , (21) 

expressing the conservation of the integrals of motion. Here the subscript "0" refers to the 
initial state. As a result of solving Eqs. (21), each parameter a, or 7^ (n = 0, 1, 2, . . .) is a 
function of the infinity of the initial invariants Eq, I^q and J„o- It is emphasized that there is 
no arbitrariness in the temperatures characterizing the Gibbs distribution of a closed system. 
In fact, this is the central problem in our theory how to determine these temperatures in 
order to predict the final state from the given initial state. 

The angular brackets in Eqs. (21) denote the ensemble averaging, 

which involves functional integrals over the space of the system states. This kind of integrals 
do not always exist. However, when the probability functional P is Gaussian, the functional 
integrals belong to the important class of Wiener integrals (their complex counterparts 
are known as path integrals [0), which are soluble and well-behaved. This is exactly what 



we use in order to resolve the ultraviolet catastrophe. In fact, we seek the long evolved 
state in the form of a coherent structure plus small-amplitude fiuctuations. This allows us 
to expand the integrals in the exponential (20) about the coherent structure up to quadratic 
terms, which will result in a Gaussian probability functional. 

The uniform and additive (with respect to the eigenmodes) invariants is another as- 
sumption lying behind the transition from the microcanonical (19) to the canonical (20) 
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distribution functionals. The additivity of the invariants can be achieved by the procedure 
of diagonahzation, which only works for quadratic forms. 

In the spirit of the conventional statistical mechanics we call the state specified by the 
detailed list of variables {fi, ■ ■ ■ , /n) the microstate, whereas the union of all microstates 
with the same (/i, . . . , /jv), N' <^ N the macrostate. The entropy S — a functional of 
the macrostate — is then introduced as the logarithm of the number of various microstates 
corresponding to the given macrostate. Up to an additive constant, S is the logarithm of 
the microstate phase volume on the manifold (19): 

'S'[/i, . . . , /iv'] = Iny" PMcifi, ■■■Jn] dfN'+i ■■■dfN ■ (23) 

In other words, the entropy is simply the logarithm of the canonical distribution functional 
(20), 

= In P[a,V;] = -a{E[a,ilj\+ Ip[a\ + J^[a,ilj\) , (24) 

if A^' 3> 1. In Eq. (24), the integrals and are defined by Eqs. (6) and (7), respectively, 
and the functions 

/?(a) = 53/3„a" and 7(«) = E7na" (25) 

n n 

can be regarded as "topological temperature functions." 

It is emphasized that the macrostate specified by N' degrees of freedom can be arbitrarily 
detailed, as we may let A^' — > oo (while preserving the requirement N' <S N), so that 
formally there is little difference between macrostates and microstates in a continuum system, 
although the apparatus of the canonical distribution (20) and the entropy (24) is much more 
convenient than that of the microcanonical ensemble (19). 

IV. B Coherent structure: the most probable state 

Maximizing the probability (20) or, equivalcntly, the entropy (24) yields "the most probable 
state" of the system. Upon varying S with no restriction on the field variations 6a and 5ip 
we obtain 

5{E[a, ij] + Ip[a, if)] + J^[a, z^]) = , (26) 
which results in a stationary solution (as(x), '^^(x)) satisfying 

= /3'(a,)+7'(a,)V27(a,) , (27) 
= -l{as) (28) 
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[compare with (9) — (10)]. There is nothing surprising in that the most probable state is 
stationary, because varying a hnear combination of the energy and the topological integrals 
(26) amounts to the Arnold (iso-topological) variation written with the help of Lagrange 
multipliers. The relation (28) stating that the fluid flow is along the magnetic field lines is 
characteristic of the "dynamic alignment" developing in course of turbulent MHD relaxation 



Similarly to the transformation (11) — (12), we can rewrite Eqs. (27) and (28) in the form 

V^T{as) = /3'(a,)/r'(a.) , (29) 

where 

r(a,) = r da^l - [7'(a)]2 . (30) 



This representation of the coherent structure will be used in Sec. |V|. 
The quantity 

7 a, X = -— — = -— — = 31 

is the local Mach number of the fluid flow. Although some interesting phenomena may 
occur near the lines where |7'| = 1, we will restrict our attention to the sub-Alfvenic case 
1 7' I < 1. A sound motivation for this is found in the absence of maximum-energy states in 
2D magnetohydrodynamics (see Appendix A) and the relaxation of turbulence to minimum- 
energy states where |7'| is necessarily less than one [see Eq. (44) and Appendix C]. Even 
though the initial condition is highly super- Alfvenic, | V^^ol ^ | Vao|, the necessary magnetic 
field will be generated by means of turbulent dynamo. 

In addition to equations (27) and (28), we must require that the equilibrium state (a^, tpg) 
be actually the maximum of the entropy S. This requirement, which is pursued in the next 
subsection, means that the coherent structure must be Lyapunov stable, which is natural to 
expect of a relaxed state. Indeed, the "fine-grained entropy" 5*, as defined by Eq. (24) is an 
integral of motion playing the role of a Lyapunov functional. 

IV. C Fluctuations: the Gaussian turbulence 

Now that we have identified (or, rather, assumed the presence of) the coherent structure 
{as,ips), we seek solution to the problem (1) — (2) in the form 

a(x,t) = as(x) +a(x,t) , ?/^(x, t) = ?/'s(x) + V'(x, t) , (32) 
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where the amphtude of the fluctuations f = (a, -0) is expected (and below confirmed) to be 
small in the long-time limit. With this in mind, we calculate the second variation of the 
entropy: 

^-oij [(V5)^ + {V^f + f3*a^ + 2i{a,)a(j] (fx , (33) 

where we denote 6a = a, Sip = ip, and (3* = /3"(as(x)) + ujs{'x.)^" {as{x)) . In order for the 
fluctuations to be finite, must be negative definite. The integral quadratic form on the 
right hand side (RHS) of Eq. (33) can be represented as the matrix element (f |iy|f) of the 
linear self-adjoint tensor operator, 

acting on a pair of functions f = (a, -0) . The boundary conditions are a = (tangential 
magnetic field) and ip = (tangential velocity) at the boundary of the finite domain. 

The orthonormal set of the eigenfunctions (a^, ipm) of W provides a natural representa- 
tion of the fluctuations: 

'^(;:)-^"(;:)' 

The standard definition of the orthonormality implies 

j (OmOn + 'ipmi^n)d'^^ = S^n ■ (36) 

Upon expanding the fluctuation field 

(37) 

in a series over the complete set of the eigenfunctions, the probability distribution of the 
fiuctuations is conveniently written as 

P[f] ^ exp = exp [-| Y: >^mfl] . (38) 



m 



which is a Gaussian distribution. 

In Appendix B we discuss the Liouvillianity of the variables fm and show that the averages 
over the distribution (38) are done by replacing VoDip by Ilm^/m Eq. (22). Then the 
fundamental averages are 

{fm) = , (/m/n) = (^mn/ {olX^) . (39) 
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The eigenvalues depend on the temperature parameters a, /3„, and 7„ and, through 
those, on the initial state. However, the behavior of the eigenvalues becomes universal in 
the ultraviolet (m ^ 1) limit. In Appendix C we show that the spectrum of the matrix 
operator W is similar to that of the standard scalar Schrodinger operator f/(x) — V^, whose 
quasiclassical eigenvalues are determined by the Bohr-Sommerfeld quantization rule: 

A-TTTT) 

A±^C±^, m»fc,25, (40) 

where S is the area of the domain and kg the characteristic wavenumber of the smooth part 
of the coherent structure. (In Sec. |V| we show that the coherent structure can also have 
singularities — current sheets — which are not important in this context.) In the Schrodinger 
case the constant C in Eq. (40) is unity, whereas for the operator (34) there are two branches 
of eigenmodes with 

1 - lYlmax < C- < C+ < 1 + lYImax • (41) 

As also shown in Appendix C, the eigenf unctions behave in the ultraviolet [Wentzel-Kramers- 
Brillouin (WKB)] hmit as 

^±(x) ~ ±a± (x) , m » k^S , (42) 

and the WKB wavenumber of the mth mode is 

^ Anm/S , m > A;^5 . (43) 

The maximum of entropy < 0) is equivalent to the non-negativeness of all eigenvalues 
Am of operator (34). It appears difficult to formulate the exact criterion of the positive 
definiteness of W in the general case; however, a sufficient condition can be derived by 
applying the Silvester criterion to the integrand in (33) considered as a plain quadratic form 
of fifth order expressed through the variables a, VS, and ^ip. Then the result is 

a > , lYI < 1 , (3*{1 - 7'2) > (Vy)^ . (44) 

Under these (or perhaps milder) constraints, the coherent structure is Lyapunov stable as 
realizing minimum of the conserved quantity —S/a = E + 1^ + J^. 

IV.D Partition of conserved quantities between the coherent struc- 
ture and the fluctuations 

So the long evolved state of the 2D MHD turbulence involves two constituents, namely 
the stationary, stable coherent structure (as(x), ^/'^(x)) and the fluctuations (a(x, t), '?/'(x, t)) 
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distributed according to the Gaussian law (38). The initial state's invariants are shared 
between the structure and the fluctuations, 



Eo = Es + E , Ifq = Ifs + If , Jgq = Jgs + Jg , (45) 

where the subscripts and s refer to the initial state and the coherent structure, respectively, 
and tilde to the fluctuations. The Gaussianity and the integral sharing properties (45) follow 
from our assumption of the small amplitude of the fluctuations, which we confirm below. In 
addition, we establish that Ip — 0, which bears a useful topological corollary. 
We start with the fluctuation energy 

E^(^ f{]a + i)Ld)(f^ = \Y. (fmfn) fUmdn + '0m^„)rf^X . (46) 

Using formula (39) and ehminating jm and Um with the help of Eqs. (34) and (35), Eq. (46) 
can be rewritten in terms of only and %l)m- At m ^ k^S the principal term in the 
fluctuation energy is 

The integrand in (47) is greater than (a^ + + It'I) ^^^^ than (am + V'm)/(l~ 

As the orthonormality condition (36) then implies, the sum (47) diverges with the number 

of eigenmodes A^ ^ 1 linearly: 

~ I / I max I / I max 

Equation (48) is a remnant of the equipartition of energy between the degrees of freedom 
(eigenmodes). At finite temperature the energy would diverge as A" — > oo, which constitutes 
the well-known "ultraviolet catastrophe." However, E is bounded from above by the initial 
energy Eq. Therefore the energy temperature l/a of the fluctuations should decrease with 
the number A^ of the effectively excited modes. Prom the conservation of energy [Eq. (45)] 
we infer 

CnN 



Analogously to the fluctuation energy, Ip and Jq also diverge at constant a, as A^ — > oo. 
However, the divergence of Ip is only logarithmic in A^, because the role of small scales is 
less pronounced in the magnetic topology invariants (6), which involve no derivatives of a 
and -0. [The linear divergence of E and Jq is due to the terms (Va)^ and u — — V^V (3) 
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and (7), respectively.] Similarly to (46) — (48), we have 

I^^ll F"{a,) {a') d'^ = Y.i^j F"{as)al d^x , (50) 



and, in accordance with Eqs. (36) and (40), 



11 ^ ^f^kax yj_< 1; 1'";^^^ ^ InTV . (51) 

2a tr A„ - 87ra(l - lyUax) ^ ' 

Upon substituting expression (49) into Eq. (51) we obtain 

' - 47r 1 - N - ^ ' 

Thus we conclude that 

If ^ Q as ^ oo . (53) 

That is, in the long evolved state, the invariants Ip are exclusively contained in the coherent 
structure {as,ips), which therefore inherits the exact magnetic topology of the initial state. 
On the contrary, the energy and the cross topology invariants, due to their linear divergence 
at — > oo, are shared between the coherent structure and the fluctuations. 

Analogously to conserved quantities, we can estimate the mean square norm of the fluc- 
tuations: 

J m 2tt{1 - iTlmaxj A/ 

Thus the mean square amplitude of the fluctuations, measured in the magnetic flux and 
the stream function, goes to zero in the continuum, or the long-time limit N ^ oo. It is 
emphasized that the amplitude of all, not only higher, fluctuation modes goes to zero. 

The assumption of a ^ a^, ip <^ ips was indeed necessary for the quadratic expansion of 
the probability functional P[a,ilj] near the equilibrium leading to the well-behaved Gaussian 
distribution. In fact, the small amplitude of a and ip is not sufficient to apply the Gibbs 
formalism, because the fluctuations in VS and V?/', which enter the integrals of energy (3) 
and cross topology (7), are found to be not small. Fortunately, the quadratic expansion 
of Eqs. (3) and (7) is also valid, because the energy and the cross topology invariants are 
themselves quadratic (bilinear) with respect to the derivatives Va and Vip. This fortune is not 
extended to many other systems, most notably the two dimensional Euler equation, where 
the resulting non-Gaussianity of the fluctuations makes it difficult to draw any quantitative 
conclusions based on the Gibbs ensemble (see Appendix D). 
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V Iso-topological relaxation and topological attractors 



Now the formal procedure of varying the integral of entropy (26) [which leads to the equilib- 
rium (28) — (29)] can be rendered more physical sense. Namely, the coherent structure min- 
imizes the energy subject to the conservation of all magnetic topology invariants {Ip, VF) 
and one of the cross topology invariants ( J^) of the initial condition: 



The term "iso-topological relaxation" describes a process whereby this minimization may be 
achieved. As we are now interested in the coherent (coarse-grained) part of the MHD system, 
and the fluctuations are set aside, the discussed relaxation is no longer Hamiltonian and 
resembles that occurring in dissipative systems. In fact, the usual dissipation in macroscopic 
systems also originates from the purely Hamiltonian molecular dynamics, where one is not 
interested in the microscopic degrees of freedom. Due to the seemingly dissipative nature 
of the iso-topological relaxation, the relaxed state may be considered as an attractor. The 
qualitative arguments developed in this section predict the appearance of the relaxed state 
without solving the complicated nonlinear problem of the reconstruction of the functions 7(a) 
and (3{a) from the initial state, a task which must be complete for a quantitative prediction 
and appears to be feasible only numerically. 

Examine the equation of the coherent structure (29), or its variational form (55). Intro- 
duce the modified magnetic flux function by the ansatz 



where the function r(a) is defined by Eq. (30) and is in principle known from the initial 
condition. We note that the conservation of the magnetic topology invariants Ip for the field 
a is equivalent to the conservation of those for the modified field a'. Then Eq. (29) can be 
interpreted as an equilibrium condition for the modified magnetic field a' with no fluid flow. 
In other words, the problem reduces to the incompressible, iso-topological minimization of 
the modified magnetic energy 



starting from the specified initial condition ao(x) = r(ao(x)). Upon minimizing Eq. (57) sub- 
ject to the iso-topological variation 6a' = {fi, a'} with arbitrary fi{x, y) results in { V^a', a'} = 
0, implying a functional dependence between the modified magnetic flux and the modified 
current. 



min E[a, if)] 



(55) 



/_F[a]=const, VF, and J-y{a,il>\=const 




(56) 




(57) 
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(a) (b) 

Figure 1: Asymptotic separatrix structures resulting from iso-topological relaxation. Topo- 
logically nontrivial initial state (a) leads to the formation of a relaxed state (b) with current 
sheets (shown bold). The arrows indicate the direction of the magnetic field. 

Once the relaxed state a'ooi^) is found, the coherent structure of the original magnetic 
field is recovered by inverting Eq. (56): 

a.(x)=r-i(a'^(x)) . (58) 

Then the stream function of the fluid flow in the coherent structure is given by 

V'.(x) = -7(a.(x)) . (59) 

The kind of relaxation undergone by the modified field a' will take place in an incom- 
pressible, viscous fluid with an ideal conductivity, where the viscosity damps down the fluid 
motion. It is well known that such an iso-topological relaxation may not be attainable in 



the class of smooth magnetic fields [^, |4g, ^ . In two dimensions, these are the saddle (x) 
points of the initial magnetic field that lead to singularities — current sheets — in the relaxed 
field. It is important to note that the location and the shape of the current sheet is not 
locally determined by the x point alone, but rather depends on the shape of the separatrix 
coming through the x point. The qualitative arguments of Ref. show that each initial 
magnetic separatrix, in course of the iso-topological relaxation, turns into a characteristic 
structure with a current sheet — the asymptotic separatrix structure shown in Fig. 1. The 
orientation of the current sheet is such as to lie within a "figure eight" separatrix and to 
border the outside of an "inside-out figure eight." 
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It appears that the final state of the iso-topological relaxation is uniquely determined 
by the initial state Oq. The same is true of the corresponding magnetic fluxes and ao, 
once the function 7(a) [and thereby T{a)] is known. Even without any information about 
7(a) the appearance of the relaxed state is well understood qualitatively through the above 
construct, because applying a function to a does not change the geometry of magnetic fleld 
lines. 



V.A Comparison with numerical data 

The computation of the long-time evolution of nearly ideal 2D MHD turbulence reported 
by Biskamp and Welter ||^ clearly shows current sheets terminating at Y points, which are 
characteristic of the asymptotic separatrix structures, although the reconnection due to flnite 
magnetic (hyper) diffusivity smears out the individual topology of separatrices. 

The spatial distribution of the fluctuations (a, tp) is determined by the eigenfunctions 
of operator (34). The potential of this operator involves second spatial derivatives of 
(through the term /?*). The singularities of js = — V^a^ are delta-function singularities at 
current sheets. This must lead to the localization of the wavefunctions (the fluctuations) near 
the potential wells (the current sheets) where 7"(as) V^7(as) < 0. This kind of localization 
of the microscopic turbulence near the current sheets is indeed observed in the computation 



of Ref. 0. Earlier simulations of turbulent magnetic reconnection [0, |5T| also conflrm this 
picture. 



V.B Iso-topological relaxation and magnetic reconnection 

So far we were mostly concerned with the ideal model of two-dimensional MHD, and the 
question is in order as to the evolution of a more realistic dissipative system involving flnite 
electrical resistivity and fluid viscosity. In general, this is a very difficult problem, because no 
straightforward perturbation theory can be built for small coefficients appearing in front of 
higher derivatives in the equations. We therefore restrict ourselves to the qualitative analysis 
of the role of small dissipation. 

If the dissipation is small, the system behavior clearly must resemble, up to a certain 
point, the prediction of the ideal MHD theory. The deviation of a weakly dissipative evolution 
from the ideal behavior is always a matter of time of the evolution. In order to neglect the 
effects of dissipation in MHD turbulence relaxation, not only must the resistivity 1] and the 
viscosity u be small but also the length scales should be sufficiently large. The ideal MHD 
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evolution discussed above does lead to the formation of small scale structures which trigger, 
in the long run, the strong effects of the weak dissipation. 

If the initial state is smooth, the small-scale structures do not appear at ones; it takes 
several nonlinear (eddy turnover) times for the small scales to show up. In the meantime, 
the system evolves towards, however not quite attains, the ideal statistical equilibrium. In 
fact, the principal manifestation of approaching the statistical equilibrium is the separation 
of scales into long-wavelength coherent structures and short-wavelength fluctuations. It is 
reasonable to assume that this separation of scales is not only necessary, but also sufficient 
for the statistical equilibrium to set in. Then, by the time when the initially small dissi- 
pation becomes important, the coherent part of the turbulent field is essentially built by 
the statistical mechanics of ideal MHD turbulence. The smaller the dissipation, the shorter 
scales are allowed to evolve in the Hamiltonian fashion, and therefore the closer the attained 
shape of the coherent structures to the exact predictions of the Gibbs-ensemble theory. 

The time scale r* specifying the crossover from the ideal regime to the dissipative regime 
is certainly much shorter than the diffusive time ~ {kgC^ri)~^ and may not be very long 
compared to the characteristic nonlinear time ta- Numerical results [§, Q indicating the 
enstrophy decay in 2D fluid in just a few eddy turnover times suggest that t*/ta is a small 
power or even logarithm of the large Reynolds number. The fast crossover to the dissipative 
regime directly indicates the fast production of small scales and, therefore, the equally fast 
approaching to the statistical equilibrium. 

After the approximate equilibrium is set in, the dissipation takes over and the small-scale 
fluctuations are significantly damped over several times r*, whereas the coherent structures 
remain little affected, at least in the case when these structures involve no singularities. 

If the initial magnetic field has x points, the coherent structure will develop current sheets. 
The coherent structure will then undergo fast magnetic reconnection. The reconnection 
occurs in a characteristic time much longer than the Alfvenic time ta, if the magnetic 
Reynolds number Rm = Tri/rA is large. By different models, Tr/rA ranges from Rll^ m, H 



to {InRmY, p > |5^, although the former (Sweet-Parker) model appears to be more 
typical [p^ . 

So the ideal MHD turbulence theory describes the early, t < min(r*,rr), iso-topological 
stage of the turbulent MHD relaxation and predicts the appearance of the coherent structures 
entering the later stages where magnetic reconnection and/or viscosity play the dominant 
role. Even then, some topological invariants survive better than others, also providing useful 
variational tools for the prediction of fully relaxed [0, ^ |59|, ^ or selective-decay ||, |[ 
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states. 

Our theory can be used to qualitatively describe the relatively early stage, ta <^ t ^ 
{TATf]Y^'^i of the nonlinear kink tearing mode in a tokamak, where two dimensional MHD 



models are commonly used for helically symmetric magnetic perturbations |40, Bq, 41, 32 



The kink tearing is accompanied by changes in magnetic topology. First, an x point in the 
"auxiliary magnetic field" = B — gBg is created near the linearly unstable g = 1 surface. 
Then the resulting "magnetic bubble" is pushed to the exterior of the plasma column by 
essentially ideal MHD motions. This process is likely to be of turbulent nature and, until 
very small scales are generated, the ideal turbulent relaxation will proceed in the direction 
of forming a coherent structure with a current sheet corresponding to the initial x point, 
as suggested by the Gibbs statistics. This stage of evolution may be pretty long, as the 
magnetic Reynolds number in tokamaks can be of order 10^ and more. Later on, magnetic 
reconnection via the current sheet will occur at a characteristic time of order (r^r^)^''^. 
Dynamically, the reconnection develops through a sequence of singular MHD equilibria with 
the same local helicity as analytically described by Waelbroeck |6^. The first of the 
sequence of these current-sheet equilibria can be interpreted as the asymptotic separatrix 
structure arising from the initial state via the iso-topological turbulent relaxation. 



VI Summary and conclusion 

The main result of this paper lies in working out the Gibbs statistics for a Hamiltonian PDE 
system with an infinity of constants of the motion. This formalism was demonstrated in 
the example of two dimensional magnetohydrodynamics but can be carried over to other 
systems. We review again the principal steps of our approach in terms of a general nonlinear 
Hamiltonian system describing the fields i/'(x, t) and having a finite or an infinite number 
of invariants !['«/'] = AiV']? I2['^]: ■ ■ ■■ Here it does not matter what these invariants are; one 
can think of Ji as the energy and of the rest as topological invariants. 

(a) The solution to the underlying nonlinear system is sought in the form i/'(x, t) = i/>^(x) + 
?/?(x, t), where V's(x) is yet unspecified stationary, Lyapunov stable solution (coherent 
structure). We then anticipate that the amplitude of the fluctuation field t) is 
going to be small and hence the exact integrals of motion can be expanded about the 
coherent structure V's up to quadratic terms: I = + I. 

(b) The Gibbs ensemble is introduced in the fluctuation space in the standard form of 
the exponential of a linear combination of all invariants, P[ip] = exp(— a ■ I), where 
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a — ai, a2, ■ ■ ■ are the reciprocal temperatures to be determined from the initial state. 
Having an infinity of Casimirs, which depend on an arbitrary function, is not an obsta- 
cle, because any linear combination of the Casimirs is again one of them. In order to 
have non-diverging fluctuations, we exercise our right to suitably choose the coherent 
structure ipg. Namely, ipgi^) is required to minimize the linear combination a • 
of the invariants, where the stationarity of the resulting state is ensured by the Arnold 
variational principle. Then a • I is a positive deflnite quadratic form, and the Gibbs 
distribution of the fluctuations is a Gaussian distribution. Prom now on, the standard 
Boltzmann-Gibbs statistics is applied in a straightforward way, at least for a flnite 
dimensional approximation using N eigenmode amplitudes fm satisfying the Liouville 
theorem. 

(c) The eigenmodes are introduced such as to diagonalize the Gibbs exponential, a ■ 1 = 
{0^1/2) J^rn -^mfm- Then averages can be computed in the conservation laws, Iq = 
Ig + in order to infer the equations for the temperatures. 

(d) The fluctuations' share of the invariants, = 1/(2q;iA^)9^I/9/^, when ex- 
panded in the eigenmodes, turns out to diverge as N ^ 00, unless the temperatures 
1 / am are let to zero (even then the temperature ratios remain flnite and keep useful in- 
formation). This is the "ultraviolet catastrophe." The regularization of this divergence 
requires the reciprocal temperatures to also diverge, e.g., ai{N) oc N. 

(e) If the square norm of the fluctuations ^V'l'^^ = Z^m(<^i^m)~^ diverges at constant 
temperatures slower (e.g., logarithmically) than the fastest diverging invariant (say, the 
energy Ji), then the average norm goes to zero as 00. This is the crucial point, 
which justifies the assumption of the small amplitude necessary for the Gaussianity 
of the fiuctuations in the given representation. If this condition is not fulfilled, one 
can always pick other variables involving lower order of derivatives, such as '0' = 
V^^i/?, and repeat the above steps. However, the exact Gaussianity of the fiuctuations 
requires another important property of the integrals of motion, which is independent 
of the variables used. Namely, each invariant must be not more than quadratic in the 
highest-order-derivative variables. Then the quadratic expansion will be also valid even 
for those (fastest diverging) invariants, whose fiuctuations are finite. This property 
holds for 2D MHD, but it does not for 2D Euler turbulence or Vlasov-Poisson system 
(Appendix D). 
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(f) If, in addition, there are invariants (such as magnetic topology invariants) diverging 
slower than the fastest diverging integral of motion, then the average fluctuation's share 
of those invariants vanishes as ^ cxd. The presence of such invariants simplifies the 
analysis of the coherent structure. 

We use the above steps to study the relaxation of ideal two dimensional MHD turbulence, 
where both infinite sets of topological invariants, magnetic (6) and cross (7), are incorporated. 
We show that accounting for all topological invariants leads to the prediction that the long 
evolved MHD turbulent state consist of a coherent structure (the most probable state) and 
a small-amplitude, small-scale Gaussian turbulence (the fluctuations). The fluctuations are 
small if measured in terms of the magnetic vector potential a and the flow stream function 
ip. The fluctuations in the magnetic field B and the fluid velocity v are of the same order as 
in the coherent structure. The fluctuation current j and vorticity uj are infinite in the long 
time limit. 

We find that in 2D ideal MHD turbulence the coherent structure has the same magnetic 
topology as the initial state, while energy and cross topology are shared between the coherent 
structure and the fluctuations. Therefore, for a sufficiently wide class of initial conditions 
having the same topological invariants, the final coherent state is the same, whereas the 
fluctuations, when measured by the standard norm (54), become asymptotically "invisible." 
In this sense, the coherent structures emerging from the turbulent MHD relaxation can 
be regarded as "topological attractors," even though the underlying dynamics is perfectly 
Hamiltonian. (The theorem of the absence of attractors in Hamiltonian systems is not valid 
for infinite dimensional PDE systems.) 

The presence of the fluctuations on the top of the coherent structure is conceptually 
important even though the amplitude of these fluctuations goes to zero in the long time 
limit: the fluctuations appear as the storage of the "lost" integrals of motion, if only the 
most probable state is compared with the initial state. This explains the well-known result 
(cf. [^D, that the topological invariants of the coherent vortex emerging from 2D Euler 



turbulence, are different from those of the initial state. In 2D magnetohydrodynamics the 
role of the initial topology is more important. In Appendix D we discuss the application of 
the Gibbs-ensemble formalism to the two dimensional Euler equation. 

We formulate the variational principle of iso-topological relaxation, which allows us to 
predict the shape of the coherent structure for the given initial state. We show how the 
problem of the ideal MHD relaxation with plasma flow is reduced to the viscous relaxation 
of magnetic field with no flow in the final state. The numerical results suggest that the 



25 



asymptotic separatrix structures with current sheets are indeed observed during the turbulent 
relaxation. It appears that these structures are the route to reconnection in the nonlinear 
kink tearing mode in tokamaks. 

Many problems of MHD turbulence remain, most notably the role of small dissipation. As 
discussed in Sec. [V.B| , this is the dynamics of producing small scales which determines when 
and how the dissipative processes become important. In order to study the phenomena of 
crossover from the ideal to the dissipative turbulent relaxation, the nonequilibrium dynamics 
of the ideal relaxation must be worked out. It appears that the formalism of the weak 
turbulence theory |]T3|, 0] can be appropriately suited for Eq. ( |B.1D in order to study the 
nonequilibrium statistics of 2D MHD turbulence. However, the important role of the ideal 
Gibbs turbulence for weakly dissipative systems is found in that the ideal turbulence forms 
predictable coherent structures, which enter the later, dissipative stages of the turbulent 
evolution. 

The comparison of the MHD and the Euler turbulence prompts us to distinguish between 
three kinds of advected fields. The first kind is passive field, such as the concentration of a 
dye or the temperature that do not affect the advecting velocity field. Passive fields tend to 
become spatially uniform due to turbulent diffusion. The second kind is active field, such as 
the vorticity in Euler fluid, which does affect the velocity field but whose lines or contours can 
be indefinitely stretched at no significant energy price. The active fields therefore tend to self- 
organize assuming topologically simple structures, like monopole vortices, whose topology is 
different from that of the initial state. The third kind can be referred to as "reactive field," 
such as the magnetic flux frozen into an ideally conducting fluid. Stretching of magnetic field 
lines is energetically expensive and cannot last indefinitely. The topology of the reactive field 
is therefore much more robust than that of passive or active fields, and the self-organization 
can lead to nontrivial coherent structures with singularities (current sheets). 
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APPENDICES 



A Extremal properties and stability of MHD and Eu- 
ler equilibria: Positive and negative temperatures 

Lyapunov stability of equilibria in conservative systems depends on the extremal properties 
of their invariants. The strongest form of the relevant Lyapunov theorem states that if a 
state realizes a conditional nondegenerate extremum (that is, a minimum or a maximum) 
of one of the invariants subject to the conservation of one or more of other invariants, then 
this state is stable. It is emphasized that the above is the sufficient criterion for nonlinear 
stability with respect to any sufficiently small perturbations. Originally formulated for finite 
dimensional systems, the Lyapunov theorem is extended to PDE systems on a case-by-case 



basis and depends on the functional norm used to specify what "small" means p5 |. 

The existence of such extremal, and therefore stable, states can be detected by inspecting 
various inequalities involving integrals of motion. As an example, consider 2D MHD states 
with the fixed cross helicity (8), 



Va-V^d'x=Ji. (A.l) 

Then the energy (3) is clearly bounded from below, 

E>\Ji\ , (A.2) 

which indicates the existence of a minimum-energy state subject to the conservation of the 
topological invariants. An example of such a state is any circular magnetic configuration 
a = a{r) with no fiuid fiow, ip = 0. This configuration assumes the absolute minimum of 
energy with respect to all neighboring iso-topological MHD states. Indeed, the magnetic 
force is the tension of magnetic field lines, which tend to shrink while preserving the area 
inside them. The smallest perimeter at fixed area is assumed by a circle. 

A formal proof, which can be extended to more interesting geometries, can be given in 
terms of the iso-topological energy variation (15) — (16). The first variation is vanishing at 
an equilibrium state, and the second is given by the general form 

S^Em = y'[-{/i,a}{/i,j} + (V{/i,a})']rf2x, (A.3) 

6^Ef = I [(V5^)2 - {/i, ou} - {fi, a} - {/i, a}{u, ^}] d^^ , (A.4) 

for the magnetic and the fiuid kinetic energy, respectively. 
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Upon letting a = a{r) and Fourier expanding the azimuthal-angle-periodic displacement 
function fJ^{r,6) = Aim('") exp(im6'), several integrations by parts (presuming that /i goes 
to zero sufficiently fast for both r and r oo) yield 

6'E^= £ m2/^°^^(^^y(m2|;x„(r)p + |/iar)-r/iUr)r)>0. (A.5) 

Equation (A. 4) also implies S'^Ef > for ip = 0. That is, any circular magnetic field 
without fluid flow assumes the conditional minimum of energy subject to the conservation of 
topological invariants and is therefore Lyapunov stable. Adding a small fluid velocity along 
the circular magnetic field lines will not change the stability of such an equilibrium. 

It is easy to see that there is no upper bound on the MHD energy, even though all the 
topological invariants (6)- (7) are fixed: by an incompressible convolution of the magnetic 
field lines the magnetic energy can be made arbitrarily large. 

The kinetic fluid energy can be represented in the form of "electrostatic" energy of charges 
with the density u, 

= \j ^(xi)^(x2)G'(xi,X2) (i^xi(i^X2 , (A.6) 

where the Green function G [= — (l/27r) ln|xi — X2I for an infinite domain] plays the role 
of the interaction potential. The total charge, Jq = J ud^yi, is conserved, and the charges 
are allowed to redistribute only along the magnetic field lines [in MHD: invariants (7)] or by 
means of incompressible interchanges [in Euler equation: integrals (D.3)]. 

The extremal properties of two dimensional Euler fiuid are different from those of MHD 
in that the energy has also an upper bound at fixed topological invariants. According to the 
Schwarz inequality, we have 

^^(xi)c<;2(x2) d^^id^^2 J ^^(xi, X2) d^^id^^2 = const /| , (A.7) 

where I2 = J uj'^d'^^K is the conserved fiuid enstrophy, and we have used the quadratic integra- 
bility of the Green function, which has only logarithmic singularity at xi = X2. By collecting 
the most intensive charges closest to each other — that is, forming a circular vortex with a 
monotonically decreasing vorticity of the same sign, we construct the maximum-energy state 
||33| . Such vortices play an important role in 2D Euler turbulence [0. Conversely, a mono- 
tonically increasing vorticity in a circular domain assumes the minimum of energy under 
fixed topological invariants. 

Adding even a small amount of magnetic field to the maximum-energy stable vortices will 
drive them unstable, because the negative-energy waves perturbing the vortices will dump 
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their energy into the magnetic field via magnetic dynamo and thereby grow in amphtude. 
The resulting turbulent relaxation will lead to other, minimum-energy equilibria, which may 
be non-circular due to the competition of the elastic tension of magnetic field lines with the 
repulsion of the "vorticity charges" threaded onto these lines. 

As discussed in Sec. [IV. C|, the Gibbs distribution of fluctuations about an equilibrium 



coherent structure has the form of 

exp[— Q;(5^(£' + Topological Invariants)/2] , (A. 8) 

where a = 1/Te is the reciprocal temperature conjugate to energy and 5^ the second (un- 
constrained) variation specifying the fluctuations. In order for the fluctuations to be flnite, 
the temperature must be positive for minimum-energy states, as it is in ordinary world or in 
2D MHD, and negative for maximum-energy states, as is possible in 2D ideal fluid ^ ^ 



B Liouvillianity of dynamical variables 

The probability distribution functional (38) is meant to describe the statistical properties of 
the fluctuations. The dynamics governing the amplitudes fm{t) can be written by substitut- 
ing the eigenmode decomposition (32) into Eqs. (1) and (2). Then, using the orthonormal 
properties of the eigenmodes, we obtain 

n np 

where 



Hmn = J [flml^n, Os} + OmlV's, «n} + 

iOmilpn, i^s} + i^mii's, ^^n} + ^m{jn, + ^m{js, dnlJc^^X , (B.2) 
Hmnp = J [am{lpn, Op} + UJ^i^n, UJp} + UJm{jn, Op}] (i^X , (B.3) 

the indices m, n, and p refer to the eigenmodes and s to the coherent structure. The diver- 
gence of the phase-space flow deflned by Eq. (B.l), 

Hmm ~l~ ^^(-^mmn ~l~ H^^mijfn i (■^■'^) 
m m mn 

must be zero for Liouvillian variables. This is generally not true even for the flrst term on 
the RHS of Eq. (B.4). However, a linear change to new variables can kill the zero-order term 
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in Eq. (B.4). An example of new variables is given by the eigenmode amplitudes of another 
linear operator, 

such that the fluctuation part oi E + Ip + can be formally written as 

a \ 



{I 


w 







(B.6) 



Then, in terms of the new eigenmodes, the equation of motion can be written exactly as 
Eqs. (B.l), (B.2), and (B.3); however, the equations for the eigenfunctions am and Um are 
different from those of the old operator W (34). The usefulness of the new variables is that 
the zero-order term in Eq- (B.4) identically vanishes after integrating by parts for each 
m, whereas the linear term, 



} + /5*{a;m,ama„} d^x , (B.7) 



generally persists unless there is no inhomogeneity due to the coherent structure. Neverthe- 
less, the remaining linear term in Eq. (B.4) vanishes as the amplitude goes to zero. Thus 
the new variables fm are Liouvillian only in the limit of vanishing amplitudes, which is also 
the assumption lying behind the Gaussian distribution (38). The validity of this assumption 
is confirmed in Sec. |1V.D| . 

As the change from the old to the new variables is a linear transformation, the non- 
Liouvillianity of the old variables only leads to a constant weighting factor — the Jacobian — 
in the probability (38), which does not affect the averaging. Hence we may write Hm dfm in 
the functional integrals in (22) instead of VaVip. 

C The spectrum of the eigenmodes 

The eigenmode equation for the operator (34) can be written 

(A„ - /3* + V2)a^ = -I'V'^m, (C.l) 
-y - 2V7' V - (VV)] a„. = {Xm + V^)^m. (C.2) 

We are interested in the high-mode (WKB) limit m 00, where we can replace V acting 
on ttm and iprn by ik. Then we get the quadratic equation for the eigenvalues, 

\l - (2k2 + p*)Xm + k\l - 7'2) + ^k^k ■ V7'2 + k^P* + yV^Y) = , (C.3) 
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having the solutions 



= k^(l ± 7') + Y T + 0{k-') . (C.4) 



The corresponding eigenfunctions satisfy 

The ordering of the modes with number can be figured out by substituting the zero-order 
eigenfunction ratio (C.5) into Eq. (C.l): 

( ^* - V^l = -4^a . (C.6) 



VI ±7' ) 1±7' 

If 7' were constant, Eq. (C.6) would be a standard Schrodinger equation with the Bohr- 
Sommerfeld quasiclassical energy levels = (1 ± j')4Trm/S and the wavenumber (43). In 
the case of inhomogeneous 7' we arrive at formulas (40) and (41). 

D Gibbs statistics of two dimensional Euler turbulence 

Although the Euler equation, 

dtUJ = {lP,Lu} , uj = -V^iP, (D.l) 

is a special case of the MHD system (l)-(2) at zero magnetic field, the properties of two 
dimensional Euler turbulence are drastically different from those of 2D MHD turbulence. 
This is primarily due to the different structure of the integrals of motion, which for the 
Euler fiuid are energy 

E=^J{Vijfd'^, (D.2) 

and vorticity topology invariants 

If^ J F(u;) d^x . (D.3) 

The formal Gibbs ensemble 

P[ij] = exp[-a{E + Ip)] ~ P[^,] exp (~ (i'lWli^)) (D.4) 



2 

can be approximately written in terms of the eigenmodes of the self-adjoint operator 

W = -V^ + VV(^/',)V^ , (D.5) 
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whose quasiclassical eigenvalues behave as 

Attitl 



Ak^/3V, k^~^, m»fc,^5. (D.6) 



Then the fluctuation energy, 

diverges with the number of modes logarithmically, whereas the topological invariants, 

diverge linearly with A^. In this sense, we have something like the equipartition of topol- 
ogy, rather than of energy, between the fluctuation eigenmodes. The square norm of the 
fluctuation stream function, / ■?/'^(i^x, converges. 



In accordance with the principles (a) — (f) set forth in Sec. 0, we conclude that in the 
long time limit the fluctuations in the stream function —>■ 0. The fluctuation energy (and 
hence the velocity v) also goes to zero, but a flnite share of the topological invariants gets 
into the fluctuations. 

Hence the coherent structure emerging from 2D Euler turbulence has the same energy, 
but different vorticity topology, as compared to the initial condition. This state of affairs 
makes it hard to predict the appearance of the coherent structure from only the analysis of 
the integrals of motion and requires dynamical considerations. 

The real difficulty, however, is that the highest-order-derivative variable, the vorticity u, 
enters the topological invariants (D.3) in all powers, not only quadratically. This makes the 
expansion of the topological invariants (D.3) about the coherent structure inaccurate, and 
hence the Euler turbulence non- Gaussian. As a result, the formally written Gibbs distri- 
bution (D.4) is useless for making quantitative predictions, because, as mentioned earlier, 
non-Wienerian probability functionals do not allow well-defined averages. This difficulty is 
present even if the validity of the Gibbs ensemble (D.4) is not questioned. We note, however, 
that this also remains unclear whether or not it is possible to derive the standard canoni- 
cal ensemble from the first principles set forth in the microcanonical ensemble for the case 
of non-additive integrals of motion (see also discussion in Sec. |lV.Ap . Thus the "exact" 
predictions based on the Gibbs statistics of 2D Euler turbulence ||2^, ^ ^ appear ques- 
tionable, because the result depends on the discretization used to solve functional integrals, 
or, equivalently, on the relative size of the "micro-cells" used for a combinatorial treatment. 
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One example of such a dependence is given by the discretization by means of point vortices, 
where the equihbrium state depends on the (arbitrary chosen) vortex strengths. Another 
example, dealing with continuum fields, was provided by Lynden-Bell in the context of 
the Vlasov-Poisson model of a stellar system, where the equilibrium is given by a combi- 
nation of Maxwell exponentials depending on the arbitrary partitioning of the phase space 
into micro-cells. There, the source of the ambiguity is the same as in the Euler equation: 
the highest-order-derivative variable — the star distribution function /(x, v,t) — enters the 
Casimirs Iplf] = J F{f)d:Kdv not only quadratically. 

The non-Gaussianity of 2D Euler turbulence suggests that the Boltzmann-Gibbs formula 
(D.4) is not valid for such a system. Nevertheless, this does not affect the general trend of 
splitting into a stationary coherent structure and a small-scale, small-amplitude (in appro- 
priate variables) fluctuations. The dynamics of such an evolution can be also described in 
terms of structures. If the initial state has the length scale much less than the box size, the 
coherent part of the turbulence evolves through weakly interacting nearly circular vortices, 
whose number decreases with time due to the vortex merger R, 11, |66[|. 
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